Topological defects and the superfluid transition of the s = 1 spinor condensate in two 

dimensions 
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The s = 1 spinor Bose condensate at zero temperature supports ferromagnetic and polar phases 
that combine magnetic and superfluid ordering. We analyze the topological defects of the polar 
condensate, correcting previous studies, and show that the polar condensate in two dimensions 
is unstable at any finite temperature; instead there is a nematic or paired superfluid phase with 
algebraic order in exp(2i#), where 6 is the superfluid phase, and no magnetic order. The Kosterlitz- 
Thouless transition out of this phase is driven by unbinding of half-vortices (the spin-disordered 
version of the combined spin and phase defects found by Zhou), and the anomalous universal 8T c /tt 
stiffness jump at the transition is confirmed in numerical simulations. The anomalous stiffness jump 
is a clear experimental signature of this phase and the corresponding phase transition. 
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Bose condensates of atoms with non-zero total spin, 
such as 87 Rb and 23 Na Q, have recently been the fo- 
cus of intense experimental and theoretical study. The 
hyperfme degree of freedom in these systems allows com- 
plex ordered states, distinct from those occurring in more 
conventional systems of spinless bosons, and combining 
superfluidity with different types of magnetic behavior. 
The dynamics and topological defects of such states can 
be observed either destructively or in situ Q B ■ I n par- 
ticular, the experimental study of two-dimensional con- 
densates has been of special interest recently for at least 
three reasons: the superfluid transition in two dimen- 
sions (2D) is of unconventional Kosterlitz-Thouless (KT) 
type |fj (i.e., driven by unbinding of vortex defects); two- 
dimensional superfluids have power-law correlations of 
the quantum phase, rather than true long-range order; 
and two-dimensional models are appropriate for some 
current experiments 

Topological defects and ordered phases in the s = 1 
(i.e., total spin F = 1) spinor condensate have been dis- 
cussed theoretically in many papers since the work of 
Ho and Ohmi and Machida Q . Experimentally, s = 1 
systems are realized using atoms of 23 Na, 39 K and 87 Rb 
with nuclear spin I = 3/2. In this letter, we resolve the 
nature of the topological defects in the polar or antiferro- 
magnetic phase by explicitly obtaining the order parame- 
ter manifold and its first homotopy group, then show that 
a new phase results in any polar s = 1 condensate at finite 
temperature in 2D. For the polar phase, Zhou |9( found 
the order parameter manifold (U(l) x S 2 )/Z 2 by noticing 
an Z 2 symmetry omitted in earlier work of Ho 0], but 
obtained an incorrect first homotopy group 01 j that this 
homotopy group was incorrect was shown by Makela et 
al. but these authors claimed that Zhou's order pa- 
rameter manifold was also wrong. We show that Zhou's 
order parameter manifold is correct but has the homo- 
topy groups obtained indirectly in Ref. Illl 



These topological defects are crucial because they cre- 
ate a new phase in two dimensions: we find that the 2D 
polar condensate is unstable to a paired or nematic phase 
at any finite temperature. This phase has algebraic order 
in 29 where 9 is the superfluid phase, no spin order, and 
a Kosterlitz-Thouless superfluid transition driven by un- 
conventional topological defects (half- vortices). It has a 
clear experimental signature: an anomalous superfluid 
stiffness jump, which can be observed in an optically 
trapped condensate by the approach used by Hadzibabic 
et. al 12] to measure the conventional stiffness jump. 

The Hamiltonian for the s — 1 spinor condensate is Q 



H = 



where a, b, a' , b' run from m z — + 1 to m z = — I , M is the 
mass of the bosons, U(r) is the trapping potential, and c 2 
and Co are constants that depend on the strengths of the 
singlet and triplet scattering amplitudes. The matrices 
F are the SU (2) generators in the s — 1 representation. 

The mean-field ground states for U = are homoge- 
neous and unfragmented in the thermodynamic limit and 
can be decomposed using ip — y / noCi where C is a nor- 
malized spin-1 spinor pj|[21|. For c 2 > the ground state 
is "polar" , and for c 2 < the ground state is "ferromag- 
netic" ; any polar state spinor £p and any ferromagnetic 
state spinor Qp can be obtained from simple reference 
states using a phase 9 and a rotation matrix U in the 
s = 1 representation of SU{2): 



C2 



Cp = e lti U 



Cf = e w U 



(2) 



The decomposition into phase 
is not unique. 



and rotation matrix U 



2 



Topological defects in spinor condensate phases can 
be understood from homotopy groups of the order- 
parameter manifold |l3|. The order-parameter manifold 
in the isotropic polar case has been previously studied: 
the original paper of Ho Q found a U(l) x S 2 man- 
ifold, where the U(l) (sometimes written as 5 1 ) de- 
notes the manifold of values of the superfluid phase 9, 
and the 5 2 refers to directions n of the magnetic axis. 
Later work || pointed out an additional Z2 symmetry 
(9 — > 9 + it, n — > — n) leading to the order-parameter 
manifold M — (U(l) x 5 2 )/Z 2 with defects that combine 
a half- vortex of the phase with a rotation from n to — n. 

We now derive M from basic principles, showing that 
the subsequent claim of a different manifold is incor- 
rect, but we find that the previously obtained homotopy 
groups of M are incorrect and that the actual ho- 
motopy groups are those of Ref. [llj. The key result, for 
readers interested in the nematic superfluid phase, is the 
existence of phase half-vortices. The overall symmetry 
group of an s = 1 condensate is G = 50(3)s X U(1)g, 
where the subscripts indicate spin and gauge (phase) de- 
grees of freedom. First consider the simpler ferromag- 
netic case: all states in this manifold can be obtained by 
a symmetry operation G on the reference state (1,0,0). 
This state is left invariant by the group of operations (the 
residual group) corresponding to an 50(3) rotation by an 
angle a along the chosen z axis and an overall phase shift 
of —a (a e [0, 2tt)). The residual group is thus H = U(l) 
and the order parameter manifold is G/H = 50(3). 

The manifold 50(3) for the ferromagnetic case also 
appears in the dipole-locked A phase of 3 He, and the 
topological defects are well understood Consider 
the polar case: for the polar reference state (0, 1, 0), the 
residual symmetry group H has two disconnected parts: 
rotations about z, and rotations about z followed by a ro- 
tation about x (or y) by 7r and multiplication by — 1 . The 
resulting order parameter manifold G/H can be writ- 
ten abstractly as [u| M' = (50(3) s x U(l) G )/0(2) s+G , 
since H = 0(2). However, points on M' can be repre- 
sented by a combination of a direction vector n on 5 2 and 
angle 9 with (n, 9) and (— n, 9 + it) identified, so M' = M. 

We write this identification as Z^ +s to stress that it 
links the gauge and spin symmetries: the resulting order 
parameter manifold is 



M 



m x 5 2 g 



(3) 



That 7Ti for t he p olar state must equal Z was shown indi- 
rectly before [ll|, but it was also claimed that M' ^= M; 
the form for M is now used to prove explicitly that 
7i"i(M) = Z, not Z x Z2 = 7Ti(M n ) as claimed in [Icj . 
and discuss the energetics of vortices. 

Consider a closed path in M. The starting point in 
M is covered by two points of M 2 = U(l) x S 2 : pick 
one of these points and trace out the preimage in M 2 . 
Either this path returns to the same point in M 2 or to 



the other identified point. In the first case, the path 
is identified with an element of tt\{M 2 ) = Z: it wraps 
an integer number of times around the U(l) part. In 
the second case, the path wraps a half-integer number of 
times around the U{\) part, and moves from the initial 
point on the sphere to the antipodal point. Identify these 
paths with the 2-tuple (n, e) in the first case, where n is 
an integer, and (n+^g) in the second case, where g 2 = e. 
Concatenation of paths is then equivalent to addition in 
the first component, and Z 2 multiplication in the second, 
but the second entry is redundant and tti(M) = Z. 

Note that although tt\{M) = Z, just as for usual 
superfluid vortices, the fundamental vortex in M com- 
bines both half a superfluid vortex and a magnetic in- 
version. Odd vortex number in m(M) means that the 
superfluid phase winds around the vortex by an odd 
multiple of 7T, while the magnetic axis goes from n to 
— n. Only the vortices with even vortex number sur- 
vive when the anisotropy field favors m z = 0: these 
are ordinary vortices of the m z — component. Since 
U(l) x S 2 is a covering space of M, all higher homotopy 
groups are direct products ir n (M) = 7r n ([/(l)) x Tr n (S 2 ): 
tt 2 (M) = Tts(M) = Z, so in 3D there are line defects, 
point defects, and coreless structures |l4|. 

The energies of the combined spin-superfluid defects 
described by the odd elements of tti(M) can be compared 
with those of normal superfluid vortices at T = using 
(Eq. ^| . We ignore the finite vortex core energy in favor 
of the energy of the surroundings, which is log-divergent 
in 2D. The only contribution comes from the | V^l 2 term, 
and using the explicit spinors gives, for stiffness K, 

C"J d 2 x^[(Vh) 2 + (V9) 2 ]. (4) 
The energy for the ordinary superfluid vortex is 



E 



irKr(2iT /2irr) 2 = ttK log(£/a), 



(5) 



where L is the system size and a the core size. For the 
anomalous (combined) vortex, the energy cost is only half 
as large: there are now two angular variables involved, 
but each runs over an angle tt rather than 2ir. 

A starting point to treat thermal fluctuations to study 
the physics in 2D is the nonlinear sigma model (NLcrM) 
on the order parameter manifold. For the ordinary super- 
fluid, this manifold is U(l): the manifold is flat and the 
coupling is constant under renormalization-group (RG) 
transformations, which underlies the existence of a Gaus- 
sian phase and finite-temperature KT transition 0] . The 
ferromagnetic manifold 50(3) = 5 3 /Z 2 is locally identi- 
cal to 5 3 , which is curved in all three directions and flows 
to weak coupling [l^: no order is expected at T > 0. 

The polar state manifold (5 2 x U(l))/Z 2 is more in- 
teresting: locally it is the same as 5 2 x U(l), which has 
a Kosterlitz-Thouless transition for the U(l) phase but 
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no S 2 order at finite temperature. The renormalization- 
group flows (i.e., /3-functions) to all orders are unaffected 
by the Z 2 identification, but the vortices in the super- 
fluid phase now carry half the normal quantum of vortic- 
ity. Since the NLcrM analysis neglects amplitude fluctu- 
ations, it is desirable to check that such a phase actually 
exists: a prediction of the above analysis of the isotropic 
polar phase is that the finite-temperature KT transition 
should be mediated by half-vortices of the phase. 

The algebraically ordered state below the KT tran- 
sition has algebraic correlations of e 2 * e , not e 10 , as 8 is 
only defined modulo ir. Note that this finite-temperature 
sigma-model analysis for both phase and spin is distinct 
from the zero-temperature analysis of spin with quan- 
tum fluctuations |9j, |lj| and from spin nematic phases of 
bosons in optical lattices 0]: the state we find at finite 
temperature is an algebraically ordered nematic super- 
fluid with no spin order and gapless excitations. A tran- 
sition mediated by half- vortices is also found for spinlcss 
bosons near a Feshbach resonance [l8j |. 

To picture this state, suppose that e 2lB had an expecta- 
tion value rather than algebraic order: locally each com- 
ponent ip a averages to zero, but the spin-singlet combi- 



nation ^QlpQ— 2lp+l1p- 

can be nonzero: 



i = pe 1 ™ 7^ and the total current 



2iM 



M ' 



(6) 



Such a state can be thought of as nematic, since order is 
only present in 29, or paired, since the order appears in 
a two-boson operator if ip represents individual bosons. 

A direct check on the above scenario is that the 
Kosterlitz-Thouless transition will occur when half- 
vortices unbind. The KT jump in the renormalizcd stiff- 
ness or "helicity modulus" at the transition must there- 
fore be four times larger than the conventional value: 



Pc 



4P° = ^. 

7T 



(7) 



The superfluid stiffness jump by in superfluid helium 
was observed by Bishop and Reppy 19] . Recently the KT 
transition of a spinless atomic BEC was observed [l2j |: 
the prediction that an isotropic s — 1 polar condensate 
will have a jump 4 times as large is a direct experimental 
test for the superfluid nematic phase. 

Numerical Monte Carlo simulations of the polar phase 
(Fig. 1) reveal a jump in the helicity modulus compatible 
with the enhanced value (Eq. d) in the thermodynamic 
limit, and clearly distinct from the conventional jump pg. 
An additional check is provided by turning on spin-space 
anisotropy. The complete phase diagram with spin-space 
anisotropy introduced through the term 

H' = J dvg 2 ft^+i + Vl^-i) (8) 

is shown in Fig. 2. With m z = favored (g 2 > 0), an or- 
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FIG. 1: (Top) Helicity modulus F as a function of temperature 
for a set of fixed values of the parameters, a — 0.5, ao = 
4.0, T C MF = 1.0, co = 10.0, c 2 = 3.0, g 2 = (see Eq. 10) 
for different system sizes. The two lines are 8T/n are 2T/-K. 
(Bottom) The goodness of fit R 2 to Eq. 12 for different sets 
of system sizes as a function of temperature. The asterisks 
and crosses correspond to sets with system sizes 6, 8, 10, 12, 
14, 16, 18, 20 and 10, 12, 14, 16, 18, 20, 24, 32 respectively. 
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FIG. 2: The phase diagram of the ground-state of the F = 1 
spinor condensate with uniaxial anisotropy (Eq. 3). The 
order parameter manifold is shown for each phase. 



dinary KT transition occurs, and is observed numerically. 
With m z = ±1 favored (g 2 < 0), the order-parameter 
manifold is reduced to (17(1) x U(l))/Z2, where half- 
vortices of the in-plane spin are bound to half-vortices of 
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the phase, and for weak anisotropy the anomalous jump the quadratic Zeeman coupling to an external field, 
survives. The experimental source of this anisotropy is The numerics use the Ginzburg-Landau free energy 



dr 



a(V^«) W„) + a (T - T c MF )^a + ?«^V>a + %W a >**b • F a , b ,^ b + g 2 (^LV+i + V'-i^-i) 



(-2 



r 



(9) 



in which the KT transition temperature (Tkt) is renor- 
malized below by the quartic terms. The transition 
is identified via the jump in the helicity modulus T, de- 
fined by introducing a twist </> in the phase 9 in (2): 



r = 




(10) 



0=0 



Metropolis Monte Carlo simulations were performed for 
square lattices with sizes varying from 6 x 6 to 64 x 64. 10 9 
Monte Carlo steps were used at each value of temperature 
for equilibration and calculation of the helicity modulus. 

The plot of r as a function of temperature T, for C2 > 
and <?2 = (no anisotropy) for different system sizes N 
is shown in Fig. 1. A jump in T is clearly visible for 
the larger systems. The two lines have slopes 8T/ir and 
2T/n, and the jump is near the intersection of the former 
with the r curves. To establish a KT transition and 
locate T c , the finite-size scaling of r(Y, T) is investigated. 
The same procedure was used for the KT transition in the 
XY model [2(ij |. The helicity modulus scales as N — > oo 



T(N,T C ) = r c 



1 



1 



1 



2 In N + C 



(11) 



Here (= 8T c /tt in this case) is the value of the jump 
in the helicity modulus in the thermodynamic limit and 
C is a constant. Even though Eq.^Jis strictly valid only 
for large values of N, it is found to hold even for small 
values in the case of the regular KT transition. T(N,T) 
for different values of N and T is fit to the expression 
Eq.^Jwith C as an adjustable parameter. The goodness 
of fit R 2 for different sets of system sizes N as a function 
of temperature T is plotted in Fig. 2. There is a min- 
imum in the value of R 2 at a temperature T — 0.0757, 
which is identified with Tkt- The value Tkt is robust 
against the choice of different sets of system sizes N for 
the fit. A similar procedure was carried out with a value 
of Too = 2T c /tt: no such minimum was found to exist 
in the vicinity of a corresponding normal KT transition. 
This offers strong evidence of the existence of an anoma- 
lous KT transition and the nematic superfluid phase. 

To conclude, we have obtained the order parameter 
manifold and topological defects of the polar s = 1 spinor 
condensate. Using a NLerM analysis of this order pa- 
rameter manifold together with numerical Monte-Carlo 



simulations, we have demonstrated the existence of a low 
temperature spin-disordered phase with quasi-long-range 
phase order. The KT transition out of this phase, driven 
by unbinding of half- vortex defects, can be detected ex- 
perimentally by its anomalous 8T c /ir stiffness jump. 
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